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ABSTRACT 

We consider the low energy effective action of QCD below the chiral symmetry break- 
ing scale, including, in Wilson's spirit, all operators of dimensionality less or equal 
to 6 which can be built with quark and chiral fields. The effect of the residual gluon 
interactions is contained in a number of coupling constants, whose running is stud- 
ied. The resulting model is an extension of both the chiral quark model and the 
Nambu-Jona-Lasinio one. Constraints on the coefficients of the effective lagrangian 
are derived from the requirement of chiral symmetry restoration at energies above the 
chiral symmetry breaking scale, from matching to QCD at intermediate scales, and 
by fitting some hadronic observables. In this model two types of pseudoscalar states 
(massless pions and massive IT-mesons), as well as one scalar one arise as a conse- 
quence of dynamical chiral symmetry breaking. Their masses and coupling constants 
are studied. We also predict a number of low energy structural constants. We find 
out that QCD favours a low-energy effective theory which is largely dominated by the 
simplest chiral quark model, whereas higher dimensional operators (such as those of 
the Nambu-Jona-Lasinio type) can be assumed to be small. 
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1. Introduction 

Two QCD-inspired quark models are commonly usedQ as substitutes of QCD in 
the hadronization regime [1-10]. 

The first one is the Chiral Quark Model (CQM) [l| with colored quark fields 
1i) q?] i — 1) ■ ■ ■ j N c , a = 1, 2, .., N F interacting with the colorless chiral field U (x), by 
means of the Yukawa-type operator 

Ccqm = q(i$~ M (UP L + U + P R )) q + ^/ 2 tr (dpU&W 

where jd = and the projectors Pl(r) = (1 ± 7s)/2 are used. M is a chirally 

invariant constituent mass, which is assumed to acquire a non-zero value once chiral 
symmetry is spontaneously broken. U(x) = exp (m(x) / Fq) is a SU(N F ) matrix with 
generators it = 7c a T a , a = 1,...,N F — 1. This massless lagrangian (|I]) is invariant 
under SUl(N f ) (£) SU r (N f ) symmetry transformations of chiral and quark fields 

U -> tt R Uttl, q L -> Vt L q L , q R -> Q R q R , (2) 

and hence the fields n(x) correspond to the massless Goldstone bosons. For this paper 
we take N F = 2. 

The constant F is the pion decay constant, whereas f represents a bare (as 
opposed to radiatively induced) pion decay constant, which is supposed to contain 
residual gluon contributions not accounted for by the constituent mass. This bare fo 
is considerably enhanced by the radiative quark- loop contribution to a final value Fq. 
Setting fo = is equivalent to assuming that the main infrared effect of QCD — the 
dynamical breaking of chiral symmetry — is associated solely to the appearance of 
the constituent mass Mq. On the other hand, in a more precise treatment low-energy 
gluons surely contribute in other ways and shift the free-quark loop value of the pion 
constant. This effect is encoded in the bare pion kinetic term[g]. 

The second approach is provided by the Nambu-Jona-Lasinio model (NJL)|4j]. It 
contains a chirally invariant four-fermion interaction. For massless quarks it reads 

Cnjl = tq$>q + [{qqf ~ (qi5T a q) 2 } . (3) 

This lagrangian is also invariant under SU F {2) ® SU R (2) symmetry transformations 
(@) of quark fields. It is however well known that for strong enough coupling g > g c 
dynamical breaking of chiral symmetry occurs in this model. This leads to the creation 
of a massless pion state, to the appearance of a dynamical quark mass Md (similar to 
Mo), and to the presence of a massive scalar meson state with the mass m CT ~ 2Md, 
the sigma particle. 

Both models provide, after integrating out the fermions, basically the same set of 
low energy structural constants describing interactions between (pseudo) Goldstone 
bosons, in reasonable agreement with data, specially for those operators which are 



1 A direct QCD bosonization approach also exists which allows to construct light meson 
lagrangians at low energies. 



3 



quark mass independent. However, neither of these two models is really satisfac- 
tory at intermediate energies ~ 1 GeV. The CQM does not describe scalar, pseu- 
doscalar, etc. meson states with masses ~ 1 GeV, which are certainly present in 
the phenomenology [12]. The NJL model (which can be extended to include vector 
and axial- vector fields, in addition to pions and scalars||) does not incorporate the 
heavy pseudoscalar IT eitheif] with mass mn — 1.3 GeV, thus treating this channel 
asymmetrically. Furthermore, its scalar sigma meson is too light to be identified with 
an observable scalar quarkonium state. Both these models fail to correctly reproduce 
the restoration of chiral symmetry at energies just above the chiral breaking scale. 
This restoration is characteristic of QCD even at relatively low energies thanks to the 
strong power suppression of the chirally non-invariant condensates. For instance, the 
first QCD contribution to the difference of scalar-pseudoscalar two-point correlators 
is of 0(1/ p 4 ) but even this one turns out to be negligible (see section 7). 

Several ways have been proposed |2|, |^, |], 0] to modify the above models in order to 
bring them into better agreement with phenomenology and to conform to the predic- 
tions of QCD. The quasilocal approach of || (see also |J7|) represents a generalization 
of the NJL model where four-fermion operators with derivatives are incorporated, mo- 
tivated by the soft momentum expansion of the QCD effective action. For sufficiently 
strong couplings, the new operators promote the formation of additional scalar and 
pseudoscalar states. These models allow an extension of the linear a model provided 
by the NJL model, with the pion being a broken symmetry partner of the lightest 
scalar meson just as before, and with excited pions and scalar particles coming in 
pairs. 

Clearly, to deal with the problem of which is the relevant effective action of 
QCD at intermediate energies we have to think in more general terms and use 
only those constraints which can be derived from symmetry and QCD principles 
and renormalization-group ideas. We propose here a wilsonian approach. We shall 
consider the most general action describing strong interactions once the spontaneous 
breakdown of chiral symmetry has taken place at a scale A, without committing our- 
selves to any specific mechanism or model. This effective wilsonian action should 
clearly contain the new light degrees of freedom appearing below A, namely the 
(pseudo) Goldstone bosons. In this paper we thus adopt the point of view that the 
pion is a distinguished hadron state being the (pseudo) Goldstone particle associ- 
ated to the breaking of chiral symmetry. Their couplings to fermions and to other 
mesons will therefore be the ones indicated by the non-linear sigma model. As, in the 
chiral limit, the pion is the only massless hadron state one can presume that, after 
its explicit incorporation in the low energy effective action, the remaining nonlocal 
operators are expandable in powers of momenta in all color-singlet (hadronic) chan- 
nels. We thus include all possible local operators up to dimension six (the choice of 
this upper dimension will be argued latter) which can be built with the help of the 
chiral field U and quarks which are gauge and chiral invariant. Gluons themselves 
are not explicitly included. Rather the effect of low energy gluons is contained in a 
number of parameters, which will depend on quantities like (a s GG) and similar ones. 
These parameters cannot be determined at present and they will be taken as given 

2 This meson is often identified with the radial excitation of the pion and denoted as ir' . 
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and used as boundary conditions (at scale A) of our renormalizat ion-group equations. 
However, once properly normalized they are expected to be of 0(1) on naturalness 
grounds, unless forbidden due to some QCD property (such as e.g. chiral symmetry). 

Thus the effective lagrangian of strong interactions below the breaking scale will 
be taken to be 

C = q(i^-M Q {UP L + U + P R ))q+ l -flti(d^Ud^) 

+AC {5) + AC {6) . (4) 

where we have made manifest the lowest dimensional operator induced by the break- 
ing. It is the one considered in the CQM and we just remark here that is the one of 
lowest dimension and hence expected to be the most relevant one. The remaining op- 
erators have been classified according to their dimensionality and lumped into A£( 4 ), 
A£( 5 \ etc. One of these operators is the dimension 6 operator 

(5) 

but of course there are many more. Their role will be discussed later. The scale 
normalizing the 4-fermion coupling constant is supposed to be related to the non- 
perturbative chiral symmetry breaking (CSB) scale ~ 1 GeV. Namely, we expect 
that attractive 4-quark forces are generated as a net effect of high-momentum gluon 
exchange with the natural infrared cutoff induced by the scale A. We have also ex- 
plicitly included the color factor N c to allow us to use the large number of colors 
limit. In addition we work in the leading log approximation. That is, we retain only 
those contributions which are universal and regulator independent (sometimes they 
can be constant when the logarithm is absent). Technically, this amounts to saying 
that our results are exact in the limit where there is a large separation between M 
and the chiral symmetry breaking scale A. While this is of course not quite true in 
QCD, it is by no means a bad approximation and, in any case, the error made is of 
the same order as the contributions from operators which are beyond the precision of 
our model anyway. 

Let us now advance some of our results. It will be proven that the mass spectrum 
of colorless states below A (in the large- N c approach) includes two pseudoscalar states: 
the pion, massless in the chiral limit (m q = 0), and a heavy Il-meson, and also one 
massive scalar state. In our model the Il-meson, rather than the pion, is the natural 
partner of the scalar meson. We have not done so explicitly, but the model can be 
easily extended to include the lightest vector and axial-vector QCD bound states. 
Keeping operators up to dimension six thus provides a 'one massive resonance plus 
continuum' model for the spectral functions in the large N c limit. One could easily 
include excited states by considering higher dimensional operators, but at the price 
of pushing A away from the chiral symmetry breaking scale and making the meaning 
of the bare coefficients in the wilsonian action less transparent. We thus treat excited 
states as part of the continuum and integrate them out, as they are naturally heavier 
than the chiral symmetry breaking scale. (Of course in large N c QCD we can adopt the 
opposite point of view, namely to dispose entirely of fermions and gluons and assume 



y 
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that the value of the constants of the effective lagrangian at the chiral symmetry 
breaking scale A is provided by integrating the infinite tower of narrow resonances, 
or one may take a position somewhere in between.) Due to 4-quark forces the bare 
scalar decay constant verifies F a > F n and the linear a-model equality is no more 
valid for it. However the residual chiral symmetry can be implemented in a way whose 
restoration at high energies is manifest. This residual invariance provides generalized 
a model relations between decay coupling constants for three kinds of mesons. We 
will then see that phenomenology favors a weak coupling regime for the dynamics 
associated to four-quark operators. QCD seems to suggest that they play a marginal 
(albeit still phenomenologically important) role in CSB. We also derive the familiar 
chiral structural constants and touch upon basically all relevant parameters in low 
energy phenomenology. 

The paper is structured as follows. In section 2 we formulate and motivate the 
selection rules which are applied to build the most general lagrangian. In section 3 we 
classify the operators assembled in A£W according to their role in the CSB, reduce 
their variety to a minimal set with the help of the selection rules and establish what 
we call the Extended Chiral Quark Model (ECQM). In section 4 we introduce in the 
ECQM the auxiliary scalar and pseudoscalar fields in order to describe dynamical 
effects in the large N c approach. The mass-gap equation is derived and its solutions 
are analyzed in the weak and strong coupling regimes. In section 5 the mass-spectrum 
of scalar and pseudoscalar fields is obtained and the pion decay constant is determined 
in terms of the bare parameters of the model. In section 6 external vector, axial- 
vector, scalar and pseudoscalar sources are attached to the model. By these means 
the decay coupling constant of the pion is calculated and the current quark mass m q 
is incorporated. The decay constant of the heavy II is also estimated to first order in 
m q . In section 7 a large iV c analysis of two-point correlators in QCD is performed and 
chiral symmetry restoration is analyzed on the basis of the operator product expansion 
(OPE) for high Euclidean momenta. The condition of chiral symmetry restoration 
leads to two sum rules a la Weinberg, but in the scalar channel. In section 8 the fitting 
of parameters of the ECQM is performed with the help of a number of observables 
and theoretical inputs. Finally section 9 contains our conclusions. Some technical 
details as well as some considerations which are not along the main line of the paper 
have been relegated to the appendices. 

2. Selection rules 

In this section we shall present the list of the rules allowing us to select the effective 
lagrangian operators and parameters of the Extended Chiral Quark Model. They are 
the following 

1. Locality. 

2. P and C invariance of the Lagrangian. 

3. Chiral and gauge invariance. 
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4. Chiral Perturbation Theory (ChPT) counting rules for chiral fields at low en- 
ergies, and the canonical dimension counting rules for quark fields valid for 
intermediate energies. 

5. The incorporation of all operators of dimensionality smaller or equal to six, in 
order to take into consideration the qualitatively new dynamics originating from 
the local four-quark interactions. 

6. The large N c limit. 

7. The assumption that residual gluon interactions below the chiral symmetry 
breaking (interactions whose strength must, logically, be much diminished, as 
the pion carries a substantial part of the strong force) can be included as bare 
coefficients of the effective lagrangian at scale A. These contributions will be 
deemed small at this scale and hence they will only be considered for the leading 
operators in the chiral lagrangian. 

8. The subsequent neglect of bare operators of next-to-leading order in ChPT 
which contain four or more derivatives acting on mesons and no quark fields. We 
shall assume that the bulk of the contribution to these operators will be the one 
induced radiatively after integration of the fermions. (For explicit expressions 
of the bare coefficients of the Gasser-Leutwyler lagrangian in terms of the gluon 
condensate, both for the leading and next-to-leading terms, see 0.) In any 
case the eventual inclusion of non-zero bare values for the coefficients for the 
0(p 4 ) and 0(p 6 ) chiral operators adds only to the value of the Gasser-Leutwyler 
coefficients and is of little consequence for the rest of the analysis. 

9. The large cutoff approximation in the radiatively induced bosonic operators 
with the cutoff A being the scale of the chiral symmetry breaking. It will be 
applied in a strong sense as a large log approximation in the cutoff dependence 
which consistently supports the neglect of operators of dimension higher than 
6 (see!). 

10. The neglect of color singlet higher spin current-current interactions whose de- 
scription requires explicit vector, axial- vector and higher-spin bosons. Likewise, 
we restrict ourselves to the study of isosinglet scalar and isotriplet pseudoscalar 
meson physics only. The extension to other fields is straightforward. This re- 
striction does not affect the accuracy of the model other than by shifting the 
bare value of g^- 

11. The use of the equations of motion (EOM) of the CQM (with minimal number 
of derivatives) in order to replace derivatives in operators with the mass-like 
parameters and with higher-order multi-quark and multi-pion interactions. 

12. The restoration of chiral symmetry at energies higher than A, that is, the fact 
that scalar current correlators should approach the pseudoscalar ones with rel- 
ative asymptotics not worse than in QCD. 
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13. The normalization of the resulting pion decay constant F n to the phenomeno- 
logically observed value and of the quark condensate to the typical theoretical 
input (qq) = -(200 4- 250MeV) 3 . 

14. The chiral dynamics relation 

n s (o) - n-(o) = 64L 8 M^, (6) 

where llp a (0) stands for the zero-momentum value of the pseudoscalar two-point 
correlator with subtracted pion pole term projected onto the isospin singlet 



channel, and Lg is a chiral constant of the Gasser-Leutwyler Lagrangian 13 . 



3. Complete set of operators of d < 6 in the ECQM 

The complete set of chirally symmetric operators of dimension less or equal to six 
which in principle should be included in the effective Lagrangian can be divided into 
the following subsets: (a) purely chiral operators; (b),(c),(d),(e) operators bilinear 
in the fermion fields of d = 3,4,5 and 6, respectively; (f) four-fermion interactions 
of d = 6. For the description of the above operators it is convenient to rewrite the 
derivative of the chiral field in terms of chiral currents 

l^rfdJJ, r^ = -(d,U)UK (7) 



(a) The purely chiral operators relevant for the description of low-energy pion dy- 
namics include some with two chiral currents 



tr (y) = -tr (d^Ud^ 

tr(«yun, 

tr(^(5X)^ (dj*)). 



(9) 
(10) 



Similar operators with Z M — > r M are in fact not independent. In addition we have oper- 
ators with four or six chiral currents. (See and ]TJ| for a complete list.) They give 
next-to-leading contributions to strong interaction processes in the derivative expan- 
sion, namely their contribution contains at least four powers of external momenta. 
As discussed in section 2, we shall neglect these operators assuming that at the scale 
of chiral symmetry breaking A the corresponding bare coefficients are small. The 
operators are written in a form which is specially convenient in order to use the EOM 
for the chiral field. These we derive, as usual, from the lowest dimensional effective 
lagrangian, namely the CQM in (|I]), and they read 



2M, 



ft 



'0 

M , 



(<lRU)iqLj ~ qLj(U ] q R )i - -5ij{q R Uq L - q L U j q R ) 



2 \ T > i 3 



qRUrq L - q L rU ] q R 
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A similar equation^ can obtained for the right current r M with the replacement U — > 
W; Ifi —> r^, L <-> _R. Eqs. (|TT| ) allow us to replace operators @, (|T0D by four- 
fermion ones. 

(b) This class is particularly simple. There is only one operator of d = 3 and it is a 
chirality flipping one 

q R Uq L + h.c. . (12) 

(c) The operators of d = 4 are chirality preserving 

Ql + Qr fiqRi (13) 
qL Hl + qn jf-qR. (14) 

The latter operator contributes to the g& constant of the (constituent) quark JE]. 

(d) The complete set of d = 5 operators, which are chirality flipping, can be built 
from operators with a 'normally ordered' differential structure (i.e. all the derivatives 
either act on chiral fields or to the right, on quark fields Ql,Qr)- I n particular, the 
full list of these readsQ 



q R Ud 2 q L + h.c, 


(15) 


q R U(d ft l»)q L + h.c., 


(16) 


fet/P^gz + h.c., 


(17) 


qRUlpFqL + h.c, 


(18) 


q R U a" u I ^q L + h.c, 


(19) 


q R Ua^[l^l u }q L + h.c, 


(20) 


q R Ua^(d^)q L + h.c, 


(21) 



where a^ v = [7^, 7"]. These seven operators are not independent however. Namely, 
the identity 

d,A v - dj^ + [Ip, l u ] = (22) 

allows to replace ( |2T| ) in favour of (|20| ) and one is left with 6 independent operators. 
A further reduction is provided by applying the EOM, which we again borrow from 
the lowest dimensional CQM ([I]) 

% @q L = M U + q R , 1 @q R = M Uq L . (23) 

See the appendix A for a justification of our usage of the EOM. 

In order to use the simplification brought about by the EOM most effectively it 
is convenient to rewrite all the operators in such a way that the Dirac derivative acts 
either on the left or on the right side. Two independent differential identities are 
useful for this purpose 

17/0= Ul^ + ^Ua^lM, (24) 
$U / + U / = 2Ul"dn + U(d^) + UlpV 1 . (25) 

3 The last term in the first line of (|ll]) appears for SU (2) chiral fields and is absent for U (2) fields. 
4 In h.c. terms U' and right chiral currents appear in a mirror reflected way. Therefore all the 
analysis done with left currents is valid also for the h.c. terms. 



9 



Then any two of the operators ([T7|) , (|18|) , fllTf ) can be replaced with a combination of 
the third one and the operators 

q R UJ(@q L )+h.c, (26) 
q R ?{UJ[q L )+h.c., (27) 

which in turn can be related to already considered lower dimensional operators after 
application of the EOM. Similarly, the first two operators ([15|), ( |16"D can be subject to 
EOM reduction with the help of (|TTD , (|23|) . Finally we end up with only two relevant 
operators of d = 5. Let us take the operators fll8"|) and (^) as the basic ones. 



Finally, we analyze the structure of (|20"D and discover that at low energies, when 
Z M ~ id fJi iT a T a , this operator becomes 

q R Ua^[l^, l u ]q L + h.c. ~ -2%e abc d^ a d^ h q R a^T c q L + h.c. (28) 

The quark bilinear in the above equation, being proportional to the antisymmetric 
tensor current, will be interpolated after bosonization by derivatives acting on bosonic 
fields: e abc d^ a d v § h . It will involve either two pseudoscalar fields or two scalar ones 
as we do not consider any vector or tensor fields in this paper. Therefore after 
integration of the fermion fields, only terms with four or more derivatives acting on 
scalar fields will be produced by this operator. Thus ( PD| ) is on the same footing as the 
next-to-leading set of operators containing at least four derivatives of bosonic fields, 
whose bare coefficients we agreed not to consider. In fact, on dimensional grounds 
based on our selection rules, the present contribution will be suppressed with respect 
to the purely bare value of the 0(p 4 ) operators by powers of M /A and therefore is 
subleading even with respect to the bare value of the 0(p A ) coefficients, which we 
anyway neglect. 

With these arguments only one operator, say ([18]), should be taken into account 
in the ECQM. 

(e) We move next to the d = 6 operators which can be built out of a quark bilinear 
by adding derivatives and chiral currents. This category consists of 25 (chirality 
preserving) operators with three derivatives, of course not all of them independent 

QlO£ q L + q R d*q R , (29) 

where the 'normal ordered' differential operator 0% can be any from the following 

list g 

0d 2 , (30) 

fd 2 , d 2 (f), W*)?, Fd,0, (31) 

0(J M )0", (d^, 0(S M J") e^ pXl5 ^(d»F)d\ (32) 

l M l"0, 0(gi", PM, (0 M /)Z", (33) 

/(c^), (c^P)/, ^/S", fl„d», t^ pXl ,ri v l p d\ (34) 

e^ P xl^{d v l p )l\ ^ vp x lhl n v {dn x ), (35) 

W757 M W\ /l^, l^F, l^/. (36) 



3 The operators are obtained from by replacement Z M 
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We use the conventions 75 = 270717273 and e i23 = 1- Not all of them are independent, 
of course. We have the following seven identities 

d 2 (?) =d(dpi") + [/, {dpi*)] - [i", (dp/)}, (37) 

(dpf)d" = [d,lpd"]-[lp, /\d", (38) 

{l",(dp/)}=[fi,lpl"]+{l",y,lp}}, (39) 

[F,@(i»)} = r,(d»f)}-r,KAi (40) 

e,u P xl 5 Y(d u l p )d X = -e^ pXl ^ri"d\ (41) 

e^ P xl^{d u l p )l X = -e^ pXlbl n u l P l X = Wa757^(^/ A ) (42) 



and one is left with 18 independent operators in 

Of these, ten operators are subject to obvious reduction with the help of the EOM. 
They are the ones containing 

<f)d\ fd 2 , (d^l") d, l"dp 0, Ipl" 0, (43) 

w„d"] =m,)d", w^] = mwh (44) 

[0,(dpl")\, /(dpi"), (dpi")/. (45) 



Evidently the related 10 operators in (|29"D are eliminated straightforwardly. Of the 
remaining ones, three operators, namely 

[I", (dp/)], [l„Ad", (46) 
e,u P xl 5 Yl v l p d\ (47) 

are involved in the following two identities 

2[J„, /\d" + [I", (dp/)} = [/, d 2 ] + [0, (dpi")} + 2[d, Ipd"] + ff, (dpi")], (48) 
iepu P xl 5 ri v l p d x - [I pi Ad"= \op V l"l u d, (49) 

and, therefore, by using the EOM one can discard the operators (HH) in favour of 



Finally, one is left with the following 6 operators 

1 
2 



{/, Qd"-Upl"d, (50) 



e P u P x7^n p d x , 

ep Vp xl^l v l p l x , /Ipl", lp/l", Ipl"/ (51) 

Notice that one of the operators contains a piece which is manifestly reducible by 
the use of the EOM. We have kept it to make the associated quark bilinear traceless. 
Note also that the analysis of the structure of other possible EOM-reducible operators 
with a different position of the Dirac derivative d does not yield any new constraints. 

We notice that the above structures contain at least two chiral currents. Let us 
consider, for instance, the quark bilinear containing an insertion of the first of the 
above structures. Since lp ~ idpix a T a , 

q L {/, lp}d"q L ~ (-2dpK a d u Tt a + ^dpTT^TrA q L l P d u q L . (52) 
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Since we neglect tensor fields, at low energies the fermion bilinear qi^d^qi must be 
bosonized as 

A (d»d v - -^d 2 ) S + B (d^PVP - ^d p Pd p P^ 

+C (d»Sd v S - -g^dpSdPS^j (53) 

where S, P stand here for generic scalar and pseudoscalar fields and A,B,C are 
bosonization constants. As a result all fermion bilinear d = 6 operators lead, after 
integration of the fermion fields, to 0{p 4 ) terms in the low energy effective chiral 
lagrangian, which, just as before, we need not consider as they are suppressed by 
a factor of 0(Mq/A) with respect to the bare 0{p A ) terms which we have already 
neglected. Therefore there are no operators of type (e) which are eligible to be 
included in the leading wilsonian action. 

Insofar the effective operators considered are leading in the large N c expansion. 
To complete the list of chirally invariant operators bilinear in quark fields one should 
mention also those ones built of two disconnected traces in isospin indices and ac- 
cording to the QCD-motivated combinatorial estimations [TIJ they may arise in the 
next-to-leading large- N c order as a result of chiral bosonization involving two quark 
loops. For SU{N) chiral fields^] these operators are 



q R Uq L tr(l^)+h.c, (54) 

qLlA^ {IT) - \q L fetr (l^F) + (L, l„ -> R, r M ), (55) 

W?Ltr (/X) - (L, Z„ -> R, r„), (56) 

UlM^ {IT) - (L, Z„ -> R, r M ), (57) 

q L ^q L tr {IV) - {L, R,r^), (58) 

WA^L7 M gLtr (lTl x ) - (L, Z„ -> R, r„) (59) 

= -e^to^itr ((9"Z")Z A ) - (L, Z M -> i?, r M ), (60) 

<?l fetr(Z/0 + (4^ Av>, (61) 

qLl^Ltr {U&X)) + (A Z M -> i?, r M ), (62) 

qLl^Ltr {l"{dj»)) + (L, Z M -> i?, r M ) (63) 
= gz^tr (Z„(3„0) + (A ; m - A 

= -~^(^7V)tr {IX) + {L,l, - fl,r„). (64) 



The operator fl6"0|) is reducible due to (P2"|). Three operators (|0T|) , (|6"2"D, (|5"3|), ( |6~4"D 
can be eliminated by means of the EOM. The operator (|58|) vanishes in the SU{2) 
case but may contribute for any chiral symmetry group with nontrivial symmetric 
structure constants, in particular, for U{2) and SU{3) chiral groups. An analysis 

6 If the Ua(X) extension of the ECQM were considered one would have to add 32 more independent 
operators. However in this paper we adopt the viewpoint that the U(l) pseudoscalar boson is not 
chiral as it remains massive in the chiral limit. Therefore it should arise from an effective quark 
sclfmteraction. 
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similar to the one already performed leads to the conclusion that all the operators 



(55) - ( |59| ) contain at least four powers of momenta (four derivatives) and should, to 
be consistent with our selection rules, be neglected. Thus only one operator, fl5"3|), 
should be included in the leading-order lagrangian. But for the 577(2) case it coincides 
with (|T8| ) and does not introduce any new phenomenological constant. 

(f) We finally turn to d — 6 four-fermion operators. Amongst them, according to our 
selection rules, we retain only those built from color singlet scalar and pseudoscalar 
quark densities supplemented with chiral fields. We shall choose a basis for these 
operators consisting in bilinears with a definite isospin in the color-singlet channel 

$i = qLU ] q R q R Uq Ll (65) 

$2 = q L T a tfq R q R UT a q L , (66) 

$3 = qLU ] q R q L U ] q R + q R U q L q R U q L , (67) 

$ 4 = q L r a U^q R q L T a U^q R + q R Ur a q L q R Ur a q L . (68) 

where color indices are contracted within each quark bilinear. 

Depending on the choice of the chiral group, there are a few invariant combinations 
which do not depend manifestly on chiral fields. In the particular case of SU(2) there 
are two invariants 

(qq? - {ri^qf = 2$i + 2$ 2 + $3 - $ 4 ; (69) 
{qimf - {qrqf = -2$ x - 2$ 2 + $3 - $4, (70) 

while in general only their difference is invariant and [/-independent 

■iqM>q',. = *i + $2. (71) 

As discussed in section 2, for simplicity we shall only consider the isosinglet scalar 
and isovector pseudoscalar channels. (The generalization to other channels is straight- 
forward and it does not modify in any essential way the results presented in the coming 
sections.) Then, of all the four fermion operators in (p5p - (p8|) , we need to take into 
consideration the following combinations 

$50 = {q L U ] q R + q R Uq L ) 2 = 2$j + $3, (72) 
$pi = (-q L TU^q R + q R Urq L ) 2 = $4 - 2$ 2 . (73) 

There are theoretical arguments (see appendix B) suggesting that the coefficients of 
these operators, gso an d gpi, must be identical at the chiral symmetry breaking scale 
A. However these operators may acquire different coupling constants below the chiral 
symmetry breaking scale due to the use of the EOM on other operators of d < 6. 
Indeed, the reduction of the dim-6 operator @ and of the dim-5 operator ( |16"D with 
the help of the EOM (|TT|) yields a contribution into the operator ([Tj|) changing the 
value of gpi. Thus, even if the constants gso and gpi are initially saturated by equal 
high energy gluon exchange (as chiral symmetry holds in QCD perturbation theory) 
they are split after the EOM are applied. We shall keep them as independent coupling 
constants. 
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4. Minimal ECQM and mass-gap equation 

We shall from now on work in euclidean space-time, making the appropriate Wick 
rotation back to Minkowski space when required, and regularize the effective action 
obtained after integration of the fermions by using a manifestly chiral invariant cutoff. 
Let us now put together all the operators which have survived the previous analysis. 
In the chiral limit we have 

Jo 



Cecqm = tq{pA + M (UP L + U + P R ))q-^-tr(l^ 
9s 



4 - v-m-m; 

+ ^fj^ (QlU^Qr + QnUqi f ~ {-q L rU^q R + q R Ufq L 



A 2 



qRUl^lf.qL + qLU^r^r^qR , (74) 



where 



p A =? + \6g A (JP L +jf-P R ). (75) 

As just discussed, gso and gp± need not be equal. The reason to call the coupling 
constant in the last operator 5fo is that it leads, after integration of the fermions, 
to a redefinition of / . This operator is thus inessential in the analysis of two point 
correlators in the chiral limit and becomes significant only by providing a contribution 
of 0(p A ) and beyond. This contribution is suppressed by a power of Mq/A 2 with 
respect to the bare value, which we consistently neglect anyway. Away from the 
chiral limit Sfo should be kept, however, as we will discuss in section 6. At this point 
we notice something remarkable. If we ignore all bare 0(p 4 ) contributions to the 
low energy effective action, the most general low energy action looks like a hybrid 
between the Chiral Quark and the Nambu-Jona-Lasinio models. However, as soon 
as we are interested in contributions of higher order in the momentum expansion, a 
number (not too large, fortunately) of bare couplings do contribute, modifying the 
dominant contribution from the fermion loop. 

In order to simplify further the calculations it is helpful to introduce the fields 

Ql = £qL, Q R = ?q R , Z 2 = U, (76) 
Under chiral transformations 

f — hgnt = n R £hl (77) 

In these variables the ECQM action reads 

f 2 

Cecqm = iQ {fl + M ) Q - ^tr (a^) 

+ ^j~ 2 (QlQr + QrQl? - ^5 {~Q L rQ R + Q r tQ l ) 2 
+i^Qa^a^Q, (78) 
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where 



Q = Q L + Q 



R, 



1 1 

TP =fi + - tf+ -^ 5 g A 4, g A = 1 - Sg A , 

v, = ?(d4) - {d&tf, a„ = e f (^) + (d^)?. (79) 

Let us now introduce auxiliary scalar, E, and pseudoscalar, II = Il a r a , fields in 
order to make the lagrangian (|74|) bilinear in fermion variables. We can then replace 

■ /S '° (QlQr + QrQl) 2 - {-QltQr + Q R rQ L f (80) 



4AT C A 2 4iVcA 2 
by 

iQ*(£(s) + ifi(a;))Q £ + iQ L &{x) - iti(x))Q R + N C A 2 + SJ^MT] (81) 

V gso gpi J 

and an integration over the auxiliary variables. The tildes on the fields £ and IT 
remind us that they are for the time being just auxiliary variables. The shift 

E + ifi = H(x) - M - -j^a^a^ H{x) = E + ifi, (82) 

further simplifies the dependence on chiral fields and eventually one comes to the 
following lagrangian 

Cecqm = iQ(p + H(x)P L + H + (x)P R )Q + £ bos , 

C bos = ^ti(H(x) T r 1 H(x))-£ti{aX), (83) 



where 



h(x) - m - m. 



H(x) 



H^{x)-M Q - s -^a lx a, 



with 



r 1 = ( 9 i_ 9 ± ) > (84) 

_^ = J 1_ .1 = ^ + A_ ( 85 ) 

After integrating over fermion fields one arrives to the following effective action 
for the colorless bosonic variables £, £, IT, 

w =S^- n ^{ p+h(x)Pl , +h+(x)Pk ) ■ < 86 > 
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Here the Tr operation stands for the total trace in coordinates and matrix indices. 
Notice that since /q is known to be proportional to N c , the whole effective action is 
proportional to N c . In the large N c limit the corresponding functional integral will 
be saturated by its extrema. 

Our starting wilsonian action had all modes above the chiral symmetry breaking 
scale A integrated out. The information about higher frequencies was contained in the 
coefficients of the intervening local operators. Obviously, to avoid double counting, we 
must now integrate the modes of the fermions only up to the scale A. This is denoted 
by the subindex A in the trace in (|8"6|). We shall in a few lines discuss how we do so. 
Notice also that A is also the dimensional parameter normalizing higher dimensional 
operators in the ECQM action (|71|). The scale [i is an arbitrary normalization scale for 
quantum contributions to the fermion free energy and is also an essential ingredient 
in the preparation of the low energy effective action (see below). 

The Dirac operator appearing in (^), whose determinant we must compute, is of 
the form 

D = f - 75 A) + i(S + i l5 P). (87) 

Given that we do not consider here anomalous, parity odd processes involving pseu- 
doscalar fields, we only need the real part of the determinant, which is just the square 
root of the determinant of the positive operator 

bfr = - pl + l5 {p VjA} + s 2 + P 2 -z l5 [S,P] 

+ p v S - i{A, P} - 7s {A S} + 2 7 5 [Pv, P\. (88) 

where Dy =$+ f . We can now define our regulator 



ReTr log \- j = -Tr (0 (A - DD*j log j j , (89) 

where Q(x) is the step-function. 

Note that the cyclic property of the trace operation allows us to prove the in- 
variance of (|89| ) under local chiral transformations of quark fields and/or external 
sources 

D(V, A, S, P) = ttb(V Q , A n , S n , P n )tt\ (90) 



fnbb^(v n , A n , s n , p n )n^\ (bb\v n , A n , s n , p 



Tr log = Tr log I : 

V /A V ^ 

where 




n = p L n L + p R n R n = p L n R + p R n L , 
L n = v u + A n = n t ^ + A ^ L + n[d,n L , 

R% = -a% = n R (v, - A,)n R + n^n*, 

S n + iP n = Q R {S + iP)tt L , 

S n - tP Q = Q[{S - tP)Q R . (92) 
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Note that QQ^ ^ 1. In particular, one can choose the transformation (|76|) with Qr = 
£, Q L = ^ and then the Dirac operator in ( |33"D is reproduced when substituting in 

(0) 

V„ = -v ll , A lt = ~g A a fl , S + iP = H = £ + ifi. (93) 

The determinant of l) will be expanded in inverse powers of the physical scale 
A. In addition to A we have two more scales, namely the subtraction scale \x and, 
since chiral symmetry is broken, we also expect (S) = (H) ^ 0. The scale /i can 
appear only within logarithms and is normalized either by the scale (S) acting as an 
infrared cut-off, or by the scale A acting as an ultraviolet cut-off. Within each order 
in the expansion in inverse powers of A we in turn expand in inverse powers of (S). 
It can be shownf] that for A ^> (S) and for \i ^> (S) the leading contribution of this 
expansion is given by 



N. 




J <PxC X - loop = -^Tr (e (A 2 - D&) log i 

J {2A 4 (in £ + i) - 2A 2 (in g + l) tr (S> + P^ 



16tt 

+ ln-^tr ((S 2 + P 2 ) 2 -[S,P} 2 



(S) 



+[D*, S] 2 + [D*, Pf - {A„ S} 2 - {A„ P} 2 
-2i[DV, P]{A„ S} + 2i[DV, S]{A„ P} 

-l(Kf + (F«f))}, (94) 

where we have kept, of all terms of 0(A°), only those which are logarithmically 
enhanced. We thus neglect terms of 0(p 4 ) which are not logarithmically enhanced 
and terms of 0(p 6 / (S) 2 ) or 0(p 6 /A 2 ) or smaller (neither of the latter have any logs). 
The universal part of the bosonized quark determinant is given by the logarithmic 
terms, whereas the coefficients of next-to-leading terms do substantially depend on 
the shape of the regulator. In the retained terms all the information concerning the 
preparation of the effective action is contained in only two parameters: A and \i. 
The approximation based on retaining the action ([94]) is provided by the large-A and 
large- log limit. Of course we could well have used another regulator instead of the 
G function, such as e.g. an exponential damping factor, to compute the fermion 
determinant. If so the terms not included in (|9~4]), which are non- universal, would 
change, but so would the bare coefficients of our effective action at scale A, since they 
are defined by integrating out modes above A. The result should indeed be scheme 
independent. 

After substitution of ( |9~3f ) in ( 94]) one obtains the bosonic form of the ECQM 



action (^) with the lagrangian: 

£(£,£, ft) = + A-toop. (95) 

7 Compare with |lq| where a different regularization was used. 
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The saddle point of the action fl86|) obtained from the above lagrangian is the relevant 
configuration for perturbing around in the large N c limit. Due to the chiral symmetry 
the effective action does not depend on constant chiral fields £o- The effective potential 
for constant (H) = H = S + HIqT 11 can then easily be derived 



y e// (E ,n^,M ( 



o 



f A 2 
N c \ —(So 

[gso 



+ \H Q \ 4 In 



M 

2 



a\2 



H \ 2 



A 2 

) 2 + — (ng) 
gpi 



in' 



-2A 2 |# n | 2 ln^- 



E 

A 2 



(96) 



where |flo| 2 = Eq + (ITq) 2 . In the above expression for the effective potential we have 
already retained the universal logarithmically enhanced terms only. The finite con- 
stants (except one which has been kept for later convenience) have been absorbed by 
a finite redefinition of the appropriate coupling constants. The saddle-point equations 
(or mass-gap equations) will determine the extremum of V e ff 



dV f 



eff 



<9E r 







A 2 
gso 



Mo 



4vr 2 



A 2 In 



E 

A 2 



florin 



Iflol 



(97) 



and 



dV„ 



eff 



dm 







A^ 
gpi 



An 2 



A 2 In 



A 2 



Iflnrin 



Ifln 



0. 



(91 



Evidently, the above effective potential has to be renormahzed to eliminate the 
dependence on the arbitrary scale \x. To eliminate the \x dependence we redefine the 
four fermion coupling constants and the chiral invariant mass Mq 



1 



1 



0so (A) 9sM 



4tt 2 A 2 ^ 



1 



</Pi(A) 9pM 



J-i E 

4vr 2 A 2: 



M (A) = Af (/i) 



05o(A) 
050 (//) ' 



or, in differential form, 



dg, 



so 



dfx 



050 

2tt 2 



^0 



pi 



0pi 
2tt 2 



-050 



Mo 
2tt 2 



(99) 



(100) 



(101) 



These renormalization group equations coincide with those obtained by diagrammatic 
methods in the large N c limit. They are not enough, however, to fully remove the \x 
dependence. Some pieces not originally present in the ECQM are needed to render 
V e ff truly /i independent. They are related to quark self-interaction operators of 
dimension higher than 6. In particular, the logarithmic terms proportional to \H\ 4 
in fl96|) are related to eight-fermion operators of dimension 12 which generate similar 
\H\ A - terms after bosonization (see |§). In the bosonic language this implies adding 
a piece A(/z)|if | 4 in V e ff where 

dX 1 , . 

vir = ~~a 9 ■ 102 



4tt 2 



18 



That such higher dimensional operator (in the original fermionic variables) is required 
is due to the fact that chiral and canonical dimensions do not coincide. In addition 
we need to add the field independent pieces 

J^l _iv c A 4 ln^- 103 

to the vacuum energy. Then dV e ff/dfi = 0. 

Note that the renormalization group does not uniquely fix the form of the effective 
potential. We can in fact add to the low-energy effective potential V e ff any dimension 
four polynomial constructed from the renormalization group invariants A 2 and M /gso 
and still have dV e ffjd\i = 0. These pieces would correspond to bare contributions 
which are not logarithmically enhanced. 

Matching with QCD suggests that just below the CSB scale the coupling constants 
g so, 9 'pi, A, etc. must be in some sense small (for M this means M <C A). We also 
expect the total constituent mass to satisfy (H ) <C A. We would expect these 
constants to grow as we move away from the scale A and we depart more and more 
from perturbative QCD. Indeed, the previous renormalization group analysis shows 
that these constants run and, for instance in the case of gso, gpi, and Mo grow as the 
energy decreases, which is just what we naively expect. 

After renormalizing all the coupling constants to the value they take at the scale 
A, the effective potential in the leading log approximation takes the form 

y e// (s ,ng,M ) ~ V eff (0,0, M ) 

+N C { — (S 2 - 2M S ) + — (IP,) 2 

L gso gpi 

where we have neglected the term proportional to A as well as any remaining non- 
logarithmic pieces. Unless explicitly stated otherwise, from now on we will take our 
coupling constants to be defined at scale A; that is, at the CSB scale. 

We expect the pseudoscalar coupling g P1 to be bounded in its negative values 

A 2 \H \ 2 A 2 , . 

g P1 4vr 2 |# | 2 

Indeed if the opposite inequality held then the dynamical breaking of the SU(2) 
isospin symmetry down to a U(l) subgroup would occur and (II 3 ) = IT 3 , ^ 0. As it also 
would break P we cannot accept to have this phase in the hadrodynamics at normal 
conditions. If (|105| ) holds then the real scalar v.e.v. (H) = Eq, (IP) = delivers the 



appropriate solution of the mass gap equations and determines the minimum of the 
effective potential. 

As for the effective scalar coupling it may take a priori any values. Let us first 
consider the case where g$o ~ 0(\). In terms of the effective coupling the mass-gap 
equation reads 

A 2 V 3 A 2 

^(S -M ) = -^h4 
9so 4tH 



— (Eo-MoH-^ln^. (106) 
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Just as we expect M <C A, we presume that the total constituent mass S <C A. 
Solutions to the gap equation exist for g S o > 

Eo = M fl + of^n, S <M . (107) 

For gso < this solution does not provide a minimum. On the contrary solutions for 
negative gso exist for gso large enough, which is best studied by introducing 

M 2 M 2 1 A 2 

We will call from now on strong coupling regime the one characterized by gs — gp — 
0(1) and, accordingly, the weak coupling regime will correspond to gs — gp <C 1. 
The effective potential written in these variables reads 

V eff (E ,U a ,M ) ~ V eff (0,0,M ) 

+N c d^l - 2M S ) + —(Til) 2 + ^o| 4 ). (109) 

[9s gp 2 ) 

It is also convenient to introduce the variable x = So/Mo. With its help the mass-gap 
equation, retaining the logarithmically enhanced terms only, looks as follows 

g s x 3 + (x - 1) ~ 0. (110) 

This expression shows that gs (and gp) are the natural variables in terms of which 
we can speak of weak or strong four fermion interactions. Notice that the original 
variables gso and gpi were accompanied only by the A^ c and A factors to make them 
dimensionless and of 0(1) in the large A^ c limit, but it was not clear which was the 
natural normalization. For positive gs it is evident that the equation ( |110[ ) has only 
one real, positive solution < x < 1. Its asymptotic behaviours are 

£s<l, X~(g s y 1/3 , ^5>1, (HI) 

As mentioned above, in the strong coupling regime one can reach 'overcritical' 
values of the scalar coupling, i.e. gs < 0. In this case the mass-gap eq.( |110| ) has always 
one negative solution which corresponds to the absolute minimum of the effective 
potential. The asymptotic expression for this solution is 

x~(g s )- 1/3 , £ 5 <-i; *~-|^r 1/2 , 0>£ s »-i, (112) 

Evidently the solution x < always exists whereas real positive solutions arise only for 
> gs > —4/27 and they never provide comparable minima. In all cases physically 
acceptable solutions must lie in the range 



(113) 
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We notice that the sign of the total constituent mass E and the sign (so far unspec- 
ified) of Mq are opposite in the overcritical regime in order to deliver an absolute 
minimum of the effective potential. 

5. Mass spectrum of scalar and pseudoscalar states 

In order to find the mass spectrum of bosonic states in the large N c limit we expand 
the effective action W in scalar, a = E — E , and pseudoscalar, tt, II fluctuations 
around the corresponding vacuum expectation values. 

We can read the kinetic term for the scalar meson from (51), 



where 



K s {p 2 



(]lp a(p)K s (p 2 )a(-p), 



Nr. 



(2vr) 4 
2A 2 



gso 



6E 2 +p 2 / 



The normalization of the kinetic term leads to the physical scalar field a 



a 



Its mass is 



NJ a. 
2A 2 



dsolo 



+ 6E^. 



(114) 



(115) 



lie) 



(117) 



The kinetic term for pseudoscalar states is derived by means of the expansion of 
the effective action in powers of IP and a M ~ id fJ Ti{x) / Fq. The calculation of the 
quadratic part respecting the chiral symmetry can be done directly from the action 
@), the effective potential (|96|) and the mass-gap equation ( |106| ) 



W 



(2) 



d A p 



-n {p)p — n (-p) 

-TO 

/ 2A 2 \ ~ 

+U a (p)N c p 2 I + + 2E 2 / U a (-p) 

V 9 pi J 



118) 



This form can be diagonalized and normalized with the choice of the bare pion decay 
constant 

F = V / / 2 + iV c E^i/o, (119) 
and thus one determines the physical fields tt, U 



71 



7T H ; 11 



IF 



=n°. 



(120) 



/o foVKh 

As expected one obtains one massless and one massive pseudoscalar state with a mass 

A 2 



2F 2 



lofo \gpi 



(121) 
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where the relation ( |1 19|) has been used. 

In the weak coupling regime terms quadratic in A dominate in (|117|) , (|121 



2 2A 2 2 2F 2 A 2 2 7 g so f 

m CT ~ — r ; m n ~ u r , 2 = m* — with 7 = — ; 5 = — < 1. (122) 



gsoW 11 SpiVo 2 ^ 2 <7pi Fi 







Both masses are generated at the scale A and due to Jo in the denominator lighter 
than A by a logarithmic factor. Notice that if both masses are comparable, i.e. 
m^/m^ = O(l), then /q/Fq = 0(1) and /q must also be logarithmically enhanced. 

In the strong coupling regime it is more natural to use gs and gp. The corre- 
sponding spectrum is then given by eqs. ( |117| ), ( |121| ) and in terms of the notations 
in 0, pi) 



m 2 =2M 2 (l + 3z 2 ), m^^^ + z 2 ). (123) 
9s 9s 

The ratio of the masses is controlled by the constants 5 and 7 for different values of 
gs- We notice that in terms of the rescaled constants gs, gp the major mass scale is 
M and the scale A is involved only in the definition of the parameter J . 

6. External sources, F n and current quark masses 

External vector, axial-vector, scalar and pseudoscalar sources are included in the 
basic QCD quark action with the help of the extension ( j87D of the Dirac operator. In 
the low energy effective action below the chiral symmetry breaking scale implemented 
by the ECQM, external sources are coupled to fermions in the same way, but, in 
addition, they appear through the derivatives acting on chiral fields 

d,u — d,u = d,u + [v„ u] - {^, u}, 

In = U^D^U, rv = -D^UUK (124) 

It is known that in this way the CVC and PCAC identities hold. As a consequence, 
in the rotated basis of constituent quark fields, Ql,Qr, one arrives to the following 
modification of the induced effective vector and axial-vector fields in 



A, = -\g A + + Z%£ - ^ - ?M - £4$) ■ ( 125 ) 

When substituted in (g4j) they provide the bosonized matrix elements of quark cur- 
rents by performing the appropriate functional derivatives. In particular this allows 
us to determine the weak decay constant of the pion after an axial-vector source 
is incorporated into the lagrangian. Linearizing the matrix U, this amounts to the 
replacement 

dfi — > + 2iF A^ (126) 
which produces the following piece in the lagrangian 

- *F tr (<V4) ( 127 ) 
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and thus we conclude that F n = F Q , with F given by ( |119| ). One can also check that 
the decoupling of the additional heavy pseudoscalar field II from the axial current 
is exact in the chiral limit and that the chiral symmetry is realized in the EPCAC 



way ||17|| . The identification of F n is slightly more subtle away from the chiral limit. 

Let us introduce the current quark mass m q > and, at the same time, prepare our 
effective action for the derivation of two-point correlators for scalar and pseudoscalar 
currents. These correlators can be obtained from the vacuum generating functional 
of QCD if one incorporates scalar S(x) and pseudoscalar P(x) external sources by 
functional derivation with respect to them. We supplement (|73|) with 

ACecqm = n ( m g + Z m(S + il 5 P)) <?• (128) 

The constant Z m is introduced in a chirally symmetric way to provide the countert- 
erms for the renormalization of the scalar and pseudoscalar densities (see section 7). 
The current mass can obviously be reinterpreted as a non-zero constant value for 
the source S(x) = m q + Z m S(x). Eq. ( |128| ) breaks chiral invariance, but this for- 
mally holds if we assume that the external sources rotate appropriately under chiral 
transformations 

M =m g + Z m (S + iP), M^QrMQI, q L -> Q L q L , q R ->Q R q R . (129) 

After the explicit inclusion of the Goldstone boson degrees of freedom in the ef- 
fective action below the chiral symmetry breaking scale we parametrize the residual 
gluon interactions in terms of local operators. In addition to those already consid- 
ered, there are, away from the chiral limit, quite a few more that one can write by 
including the matrix M. and using the transformations ( |1 29|) . One such operator is 
the counterpart of the term with the bare pion decay constant /o, which, as we have 
seen, plays a crucial role to make the field IT a propagating one. It is 

citr (MU^ + M^ll) . (130) 



There are several operators linear in the scalar source of chiral dimension < 4. One 
of them is 

C2tr((UiM+M*U)l^), (131) 

and it is in fact the only relevant one since tr (MW + M^U) tr (y M ) is, in fact, not 

independent for SU(2) and tr lAU^M. + A4^U)D^l^ can be reduced with the help of 
the EOM. 

There are several operators quadratic in the current masses. The relevant one is 

c 3 tr [MU^MU^ + UM^UM ] ) , (132) 
since other operators quadratic in m q are either subleading in N c or ^/-independent 

tr (MW) tr (MW) + tr (M*U) tr (M^u) , 
tr (MU ] ) tr (jtfu) , 

tr (MM*) . (133) 
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We shall not go beyond this order and, in fact, we shall consider two-point correlators 
in the chiral limit only. 

In the fermion sector there are several operators with the required invariance 
properties and linear in A4. The possible dimension 4 operators without derivatives 
are 

qRMq L + q L M j q R , 
q R UM ] Uq L + q L U^MU^q R , 
q R Uq L tr (MlP) + q L U^q R tr (UM ] 



q R Uq L tr (UM ] ) + q L U*q R ti (MU ] ) . (134) 

Those containing the additional trace will be of subleading order in the large N c ex- 
pansion and thus not considered here. Operators containing twice the source A4 give 
eventually contributions which are down by a power of £o/A 2 and will be neglected 
accordingly. 

Higher-dimensional operators with derivatives containing Ai can be reduced with 
the help of the EOM either to the operators ( |134 ) or to a set of operators which 
contribute to the next-to-leading order of the large N c expansion and therefore can 
be neglected. 

In the rotated quark basis, the only two relevant operators £1X6 cLS follows 

QrZM^Ql + Ql?M?Q r . (135) 

Their existence is due to the nonlinear realization of the chiral symmetry on the fields 
£> ^ (0) resulting in a hidden vector symmetry representation on products 

£M*£ — ► h^M^h\, 

£W — h£M£% (136) 

In principle these terms have different coupling constants unless a hidden left-right 
symmetry holds, Ql <-> Q R - In the QCD low energy effective action one would not a 
priori expect to have this symmetry. However its appearance could be revealed from 
the meson phenomenology (see section 8). 

Thus we supplement the ECQM lagrangian with the following piece 



'2 

-u 



ACecqm = i(l + e)(Q R ?M?Q L + QLZM^Q R 
\ - e) (Qr^M^Ql + Q L t*Mt*Q R 



2 

+ Cl tr (MW + M*U 



+c 2 ti[{U ] M+M ] U)l^l, 
+c 3 tr (MWMW + UM ] UM ] ) , (137) 

where e and c\,C2, C3 are real coupling constants. The above mentioned hypothetical 
left-right symmetry would correspond to e = 0. We notice that the choice of the 
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couplings is made so that the mass term corresponding to M. = m q in the lowest 
order chiral interaction, {7 = 1, coincides with the QCD mass term for low energy 
quark fields. 

After the introduction of auxiliary scalar and pseudoscalar fields just as in flSlf ) 
we perform a shift of variables to simplify the structure of the quark determinant as 
in (|82|) . Thus we define 

£ + iU = a + iU' + So - M - ((- + e)^MC + (~ - e)£M f f ) - (138) 
with a = £ — S . 

The shift ( |138| ) allows us to move the dependence on scalar sources to the purely 
bosonic part of the effective action Cbos, which is now modified by the following terms 
(we keep only those which will be relevant in the following discussion) 



V 25(50 J \ 

+ U + ^) tr ((tfM+MtylJ 



\ 25(50 . 

-^atr (MU^ + M^U) - ^atr (1,1,) 
2^50 v ' gso 

+*e^-tr (H'^M^ - £M^)) 

gpi v 7 

+ L + ^M! _ 4e^A» \ tr ^[/t^f/t + A^t^t^ . (139) 
V Sfi'so 85-pi J v / 

We notice that the difference between the two kinds of mass terms appears only in 
the pseudoscalar channel. We also notice that there are two different contact terms 
contributing to the scalar- axial- axial three point function, one proportional to C2, the 
other proportional to 5fo. In fact, as it will be discussed in the coming section in 
more detail, QCD dictates that the sum of the two pieces just above the scale A must 
be zero. The only sizeable contribution to the scalar-axial-axial vertex is the one 
induced by the exchange of the a particle, which indeed vanishes as 1/p 2 . Finally, 
the scalar and pseudoscalar correlators have also some contact terms, which can be 
read from the last line in ( |139| ). QCD has also something to say on that particular 



combination of coefficients. 

Let us analyze the influence of current masses of light quarks on the mass-gap 
solution and on the effective potential. We neglect the isospin breaking effects (with 
respect to SU(2)y) and take Ai = m q I. From ( |139| ) we immediately read the modifi- 
cation in part of the effective potential linear in m q (taking into account that (a) = 0) 

/ 2A C A 2 (S - M )\ T j , s 

AV eff = 4 Cl c —^ °J. m u q , (140) 

V gso J 

where ^ 

m q = m g cos0, cos0 = -tr (ll + . (141) 
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We shall assume, as is commonly done, that (cos0) = 1. Other choices of the chiral 
vacuum lead to different phases for M , the chiral condensate, etc. 

Just as the other parameters in the theory, we have to renormalize m q to absorb 
the logs associated to the running between the scales \i and A. Its renormalization 
can, for instance, be obtained diagrammatically and is identical to that of M 

Another way of seeing this is by realizing that to make the effective potential renor- 
malization group invariant in the presence of the soft breaking induced by the mass 
term we have to add a piece similar to the first term in ( |103| ), namely 

-^(M + r<) 2 (143) 
9so 

which, unlike (|103|) , is now field dependent. Its origin is obviously nothing but the 
terms ( |130| ) and ( |132j ). We thus see that it is necessary to include the operator ( |130| ) 
from the beginning and that 

c 1 (A) = -^-^ + g 1 (A), (144) 

where ci(A) contains other contributions which cannot be deduced from renormal- 
ization group arguments alone (and which for this reason are not expected to be 
logarithmically enhanced either). In any case the combination 

2 Cl {y)gsM + N c A 2 M (fx) = 2c l {A)g S0 {A) + iV c A 2 M (A) (145) 

is a renormalization group invariant. 

The modified mass-gap equation is (from now on we take (cos0) = 1 and thus 
= m q ) 

A 2 S 3 A 2 

— (So - Mo - m q ) = In -5. (146) 

9S0 47T Z^o 

Let us use the fact that m q is a small parameter and let us calculate the solution as 
a perturbation from the chirally symmetric one 

SoK) = So(0) + m q A + 0(m 2 ) A = < 1. (147) 

In the weak coupling regime the solution is So — Mq + m q with Mq ^> m q by 
assumption. In the strong coupling regime it is helpful again to invoke the notation 
of ( |108|) to present the mass-gap equation in the form of (|110|) 

g s x 3 + (x - 1 - -£) ~ 0, (148) 

with the approximate solution 

m q 1 m q _ S (0) 

X = X() H • = Xr> H k • -, Xn = . 149 

3-2xo Mo 3g S x 2 + l M ' M 1 } 
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This expression is valid for both subcritical, g$ > 0, and overcritical, gs < 0, regimes. 
As we shall see later on (see (|189[) ) the correct sign for the chiral constant L 5 for all 
values of gs is obtained only if the bare constituent mass M has the same sign as 
m q (positive in our conventions). Due to the mass-gap equation the total constituent 
mass S which corresponds to the absolute minimum for the effective potential takes 
the opposite sign in the overcritical strong coupling regime. Thus So is positive in 
the weak coupling regime and in the subcritical strong coupling one, but negative for 
overcritical values of gs- (See the discussion at the end of the previous section.) 

Now we derive the correction to the value of the effective potential ( PS] ) at the 
minimum due to the presence of the current quark mass ( |140|) . Substituting (|146| ) 
into the effective potential (including ( 14(J| )) and keeping terms linear in m q , one finds 
for the vacuum energy 

Ve jn = V ™ n (m q = 0) + m q U Cl + ^Eg In . (150) 

The physical quark condensate in Minkowski space, renormalized at scale A will be 

°* = '<*>- = 2^t- (151) 



Namely 



2' 



C q (A)= U^ + ^l^j. (152) 

Expanding the chiral field U around the identity and using ( |140| ) one finds the 
Gell-Mann-Oakes-Renner formula |13[ for the pion mass: tu^Fq ~ —2m q C q . The 



pion mass has the correct physical sign provided that C q < 0. Since we will see later 
that in the weak coupling regime M and, consequently, S must be positive, the 
second term in the first line of ( |152|) is positive and we conclude that C\ must be 



negative and large enough, lest the pion be a tachyon. Thus the bare coefficient C\ 
provides a contribution to the condensate which is bigger than the one induced by 
the fermion loop after integration. Taking a look at ( |143| ) we will see that in the large 



log approximation this is the case. In the overcritical regime, gso < 0, the constant 
Ci is less crucial as the accompanying term becomes negative. 

Let us now calculate the effects due to the soft CSB induced by the current quark 
mass. They arise from the perturbed solution S (m (? ) of the mass gap equation (|147|) 
in the first place, but they also modify the pseudoscalar mass matrix from ( |139| ) 



1 / 8/V A 2 

A4 2) = - hm q B (n a ) 2 + — ^-em q n a W a ) , ( 103) 



with B = \C q \/Fq. We notice that the parameter e cannot be in general renormal- 
ization independent if gp\ ^ gso- Its renormalization should be given by 



(154) 
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After the diagonalization of the kinetic term with the help of variables ( |120|) (with 
II replaced by II') the mass matrix in the pseudoscalar sector Mp is still non-diagonal 
away from the chiral limit 

M r=(r l i 2^ 2), (155) 

where the following notations are employed 
ml = 2m q B , 



. 2VW c A 2 e v^o - fi 

dl = — ; = + 



gpiBofovTo fo 



A ( 2 ln J, fo\ , 2F 2 S2\ 4A c A 2 E ^e , N c E 2 g 2 A I { 



Then, rediagonalization of (|155|) induces the coupling of the physical heavy II meson 
to the axial current from ( |127| ), something which is forbidden in the chiral limit. At 
first order in the current quark mass the corresponding decay constant is given by 

Fn = F<A s& (157) 

The mass corrections to F n will be discussed at the end of next section. 

For completeness we also show the mass shift for heavy scalar and pseudoscalar 
mesons 

ml(m q ) = m 2 a (0) + m 2 ^^ + O(mJ), 

-do 

m 2 u {m q ) =m 2 u {0) + d 2 ml + O{m 4 n ). (158) 
For SU(2) these corrections are very small. 

7. Correlators and chiral symmetry restoration 

In large N c QCD the two-point correlators of scalar and pseudoscalar quark den- 
sities are given by a regulated sum over the appropriate meson poles 



nsG° 2 ) = —J d 4 x exp(ipx) (T (qq(x) qq(0))) 



E -rrhr- + Co + c?p 2 , 

n P + m S,n 



Il¥(p 2 ) = f d A x exp(ipx) (T(q l5 r a q(x) q l5 r b q(0))) 



WE^^ + ^ + ^vU (159) 

\ n P 1 rn P,n / 



Cq and C\ are regularization-dependent contact terms. 
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The high-energy asymptotics is controlled by perturbation theory and the opera- 
tor product expansion due to the asymptotic freedom of QCD. In QCD the conven- 
tional derivative, id^ — > iD^ = id^ + G^, contains gluon fields = g s \ a G a ^ where 
tr(X a X b ) = 2S ab . In the chiral limit (m q = 0) the scalar and pseudoscalar correlators 
have the following power-like behavior |18], |19[ at large p 2 

1 + — 




8A c vr 




3p 



^(qa^qqa^q) + 0(1) (160) 



for Nf = 2, in euclidean notations. In the planar (large-iV c ) limit^ 



(n P ( P 2 ) - n 5 ( P 2 )) p2 




A S p ^ 24( ^ 2 l) ™ s Cl ~ 245 2 F 4 , (161) 



where the vacuum dominance hypothesis [jig ] has been applied and we have taken |L2 
for simplicity a s (1.2GeV) ~ 1/3. 

This rapidly decreasing asymptotics is of course a consequence of the manifest 
chiral invariance of massless QCD. When comparing (|159|) and (p.6U|) one concludes 



that an infinite series of resonances with the same quantum numbers should exist in 
every channel in order to reproduce the perturbative asymptotics. 

At scales below A, where chiral symmetry is dynamically broken, the previous 
equality needs not to be valid. However, the fact that chiral symmetry is restored 
so quickly at high energies suggests itself as a further constraint to impose on the 
ECQM. We shall demand that at the scale A, and within the precision of our approach, 
the relation ( |161| ) is approximately fulfilled. This we will call the chiral symmetry 
restoration (CSR) constraint |22| (similar CSR constraints in the vector-axial vector 
channels for specific models have been considered in |2B|, . In this section we shall 
work out its consequences in the chiral limit (m q = 0) and we will restrict ourselves 
mostly to the scalar and pseudoscalar two-point functions, although some three point 
functions will also be discussed. 

In the ECQM the two-point correlators in the scalar and pseudoscalar channels 
are obtained by performing functional derivatives of the vacuum generating functional 
Z(S, P) of the ECQM 



IW) = / *x exp(w) log Z(S, P) U P=0 
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We should mention some discrepancies in the literature flq, M concerning the signs of 



the last term in (160) and acknowledge S. Peris for bringing to our attention the paper |20J where 
the signs are proven to be correct. The correctness of the above signs can be verified by realizing 
that they provide the correct Weingarten inequality pT| between pseudoscalar and scalar two-point 
correlators when applied to the leading (logarithmic) singularity at short distances. 
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Zt +C s (p 2 ) 



p2 _|_ m 2 

n pV) = / *x e*p(ipx) J a{x) J b{Q) logZ(S,P)\s, P=0 

= * ab (^ + ^ IL ^ + C P (p 2 )), (162) 
\p z p z + m| / 

where in the chiral limit the masses mj are given by eqs. ( |117| ), ( [L21| ) and the pa- 
rameters will be derived below. Since we are well below the heavier poles (located 
above A), the functions C s (p 2 ), C p (p 2 ) will be, by construction, finite polynomials in 
momenta in the ECQM 

C V) = C S Q + p 2 C? + ... C p (p 2 ) = C P + p 2 C[ + ... (163) 

With this structure in mind the CSR constraint will basically involve the composite 
states 7r, II, a of the model and decouple from the effects of heavier degrees of freedom. 
This decoupling is always possible with a suitable choice of the coupling constants of 
local operators quadratic in the external sources. 

When imposed on two-point correlators of the ECQM (in the chiral limit) the 
CSR condition leads to 

= Of, (164) 
Z a = Z w + Z u , Z (j > Z n , (165) 
Z a ml = Z n m 2 n + A s P . (166) 

As for the constants Cf and C[, the CSR constraint is automatically fulfilled in the 
minimal ECQM after imposing our selection rules since there are no contact terms 
with two derivatives acting on M. of chiral dimension < 4. There is one operator 
of chiral dimension 4 that, on dimensional grounds, could contribute to the first 
condition (|164|) , namely tr(jAAi^, but it does not since scalar and pseudoscalar 
sources appear on equal footing in it. However, ( |164j ) does provide some information 
on the ECQM, namely 

N C A 2 Ae 2 N c A 2 . . 

c 3 + — — = 0. (167) 

8^50 8<7pi 

The constraints ( |165| ) and ( [L66|) allow to determine two constants Z± in terms of a 
third one and meson masses. In particular 

Z^m\ - A SP Z n ml - A SP , . 

£o — 5 9 — ' z n — o 9 — ■ l it:, °J 

mfi — mfi — 

PC AC dictates that = 4B 2 Fq. We can use this relation in (|168|) to find the other 
Z's. From phenomenology [12| mn — 1.3 GeV and therefore Z^mfj > Asp when 



taking into account ( p.61| ). Inasmuch as Z G is positive one derives from ( p.68| ) that 



rnfj > m 2 . 



In order to use the relations (|168|) we must find the dependence of the constants 
Zi on parameters of the ECQM. For this purpose we expand the appropriate part of 
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the ECQM lagrangian — extended with scalar and pseudoscalar sources — in powers 
of meson fields 



AC 



ECQM 



2N r A 2 



-aS 



4N r A 2 e 



pja pa 



4 Cl 



9 so 9 pi 

iV c A 2 (S -Mo) 



^a TDG 

7r r 



(169) 



where Z m was introduced in (|128|) . 

It is evident that the constant Z m is required to renormalize the operators in 
and obeys the following renormalization group equation 

Z m (n) Z m (A) 



gso(fJ') gsoW 



(170) 



We shall take Z m (A) = 1. This amounts to fixing the normalization of the external 
sources at a given scale and we do not lose any generality by doing so. This convention 
simplifies the relation between Z n and the quark condensate. Let us now diagonalize 
the kinetic term with the physical fields (|116|) , ( |120| ) and adopt the following notations 



a = 



2N C A 2 



2^NJ M 2 



gsoVNJo 9 s 

4 / jV c A 2 (So-M ) 
F a { 2g, 



2C n 



Cl 



BoF g s 



so 



F a VKToMi' 



(171) 



Then the operators linear in fields and external sources relevant for the derivation of 
two-point correlators can be written as 



AC 



ECQM 



-aaS - a(3TT a P a 



a (2e 7 + pVT^' 



U apa 



T\ Za vS - JZ n n a P a =f Jz n n a p a . 



(172) 



Notice that only the sign of the second term is fixed. The sign of the other two terms 
depends on that of g~s- From the first relation in (|168|) and comparing with ( |172| ) one 
finds that 



P = sign (g f 



1 - 



A 



sp 6F 2 



Asp 



Z^mfj 



0.03. 



(173) 



Thus, from phenomenology |12 \/3\ < 1. The relation ( 165 ) gives 

2e 7 = -pVT^lp ± sJl - (3 2 , 



(174) 

and, since \S\ < 1, it follows that \2ej\ < 1. 

We can also obtain F a from (|172j ). The scalar decay constant is defined through 
the relation (recall the equivalent one for F n ) 



2B F a 



(175) 
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Therefore 

F. = = m 

£>o9s P 

We see that the CSR constraints represent a powerful tool to reduce the number 
of input parameters of the ECQM. 

Another phenomenological condition enabling us to restrict the parameters of the 
ECQM is provided by the experimental determination of the chiral constant L 8 , which 
parametrizes the following operator in the low-energy chiral lagrangian (in euclidean 
space) 

-AL 8 B 2 tr [M ] UM ] U + U ] MU ] M) . (177) 
For U = I the above operator reduces to 

-16L 8 5 2 (S 2 -(P a ) 2 ) (178) 

In the ECQM, after imposing the condition (|167|) the only origin of this term is the 
exchange of heavy meson states. This leads precisely to the sum rule displayed in 
section 2. When combining it with the solution ( |168| ) of the CSR constraints ( |165| ). 
O and the definition (0721) for Z w one finds the remarkable relation 





(z a 






f 1 . 


1 

+ — 


64B 2 


[ml 


mfjj 


16 


\m 2 a 





A 



SP 



(179) 



which represents one of the predictions of the ECQM and modifies substantially the 
estimates of P5| . The term proportional to A$p is very small and, in practice, 
negligible 

A SP 6F 2 



Z w (ml + mfj) ml + m\ 

Thus one expects that 



< 0.03. (180) 



£s > ( 181 ) 
8mn 

This is a consequence of the CSR inequality ( 166 ), m 2 < mfj, and the smallness of 
We finally turn to the chiral constant L5, which goes with the operator 

2 B tr ((DJjyDJJ(M*U + U ] M)) . (182) 



First of all, a contribution to this low energy constant seems to come from the operator 
proportional to the coefficient (02 + j i n the lagrangian ( 139 ). For U = I this 
operator reduces to the term 

tr(£I M A M ), (183) 

which jeopardizes the correct chiral high energy behavior of the following three-point 
correlator 

n S AA = J d 4 xd A yexp(ipx + iqy) (T (qq(0)q'j 5 'j fM q(x)q'j 5 'j fl q(y))) . (184) 
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For large values of momenta one expects from the asymptotic freedom of QCD that 
perturbation theory describes the leading asymptotics. But in perturbation theory it 
is evident that the correlator ( |184| ) identically vanishes in the chiral limit as any Feyn- 
man diagrams necessarily contain an odd number of 7 matrices. It is a consequence 
of the chiral invariance under the transformation 

<?0) -> 75?(a0, q( x ) -9(^)75, S(x), P(x) -> -8(x), -P{x), (185) 

which in the chiral limit holds for the generating functional of perturbative QCD. 
Therefore the correlator (184) should vanish at large momenta and cannot have 



a constant asymptotic behaviour such as the one provided by the contact term (|183[) 



in the generating functional. We assume this behavior to hold at scale A for the 
correlators of the ECQM and impose one more CSR constraint 

c 2 + — = 0. (186) 

z 9so 

Consequently, the constant L 5 is induced solely by a exchange, both from the terms 
in ( 139 ) and from the piece linear in a in (|94D 



-^N c g 2 A EaIovtr(W. (187) 

By coupling an external axial source and repeating the analysis done in the massless 
case, we determine the shift in the pion decay constant 

'-^^-f+^W (i88) 



in the notations of ( 147 ). This implies 



U = j^ff? - fl + ^A). (189) 

Notice the appearance of Sfo in these formulae. As we announced, this effective 
constant is relevant in the m q 7^ case. This is a nuisance because it introduces a 
further unknown parameter into the problem. However, it is clear from ( |119| ) that 
Fq — /o is logarithmically enhanced, something which is not expected to happen for 
5f . In the large log approximation it is thus justified to neglect 5f in front of the 
other quantities and we will do so. 

Notice also that ( |182| ) is linear in m q and the sign of L 5 is correlated with that of 
A and E . Indeed, for any value of g$, the product g~s x is positive, while the sign of 
A is always the same as the one of x. Thus the sign of L 5 is finally that of M and 
we conclude that M must be positive, just like m q . Had we chosen the opposite sign 
for m q the sign of M would have came out reversed as well. 

8. Fit and discussion 

Let us enumerate all the parameters of the ECQM. The first seven parameters are 
not directly observable: M , x (instead of S ), 5 (instead of / ), gs, 7 (instead of gp), 
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Iq (instead of A), e. Three other parameters are directly taken from phenomenology: 
g A ~ 0.8 m, m g (lGeV) ~ 7 MeV and B ~ 1.4 GeV ||, (S„ instead of c x ). The 



two last ones provide the pion mass ~ ^2m g B ~ 140 MeV. One more parameter, 
Asp is related to the four-quark condensate which can be estimated as in (|161|) by 
employing the vacuum dominance hypothesis. However this parameter makes almost 
no influence on the main meson characteristics (as it follows from ( 173 ), ( |1 79|) the 



deviation is of order 2%). Therefore we neglect it in the forthcoming fit. 

The above 7 unknown parameters are involved in the determination of 7 physical 
quantities: m a , mn, F , F n , F a , L 5 , L 8 , but they are connected by 3 relations; namely, 
the mass-gap equation ( |11U| ) and the two CSR constraints ( |174| ) and ( |173| ). Thus 
the set of physical entries is sufficient in principle to find all the ECQM constants 
from phenomenology and make some non-trivial predictions. In a way this is a little 
miracle since we started with O(50) operators, each one with an unknown constant. 
In addition our model reproduces the same value for the chiral constants L\, L2, L3, Lg 
and L10 as the simplest chiral quark model in what concerns the radiatively induced 
quark contribution. These cannot be directly read from (p4|) because we have retained 
there only the logarithmically enhanced terms. It is known, however, that the pieces 
in the effective action leading to the above Lj's, albeit not logarithmically enhanced, 
are universal and independent of the regulator. Of course the direct bare contribution 
to the Lj's is not included in our approach (see appendix C). 

Let us perform a tentative fit using the phenomenological values for the heavy 
pion mass mn ~ 1.3 GeV [12] and for Fq ~ 90 MeV [ 13fl . First we relate the scalar 



meson mass to the chiral constant Lg by means of (|179 ) 



ml = mU 16F0 L _ i) _ 1 ^ (19Q) 

When taking the above input for mn, Fq and the phenomenological value of L$ = 
(0.9 ± 0.4) x 10 -3 [|13j one obtains the mean value for the a meson mass m a ~ 0.9 



GeV which is well compatible with the experimental data [12]. The central value 
corresponds to |/3| ~ 0.7. 



Next let us use the definition of the pion decay constant ( 119 ) in order to eliminate 
the parameter I 

FoVT^ 



NJo = 1 ■ (191) 

9a\x\M 



After substituting in (|171|) one obtains 

, oi B g A g s x 



(192) 



This relation can be combined with the definition of the chiral constant L 5 ( 189|) and 
with eq.( p.23| ) for the scalar meson mass reexpressed in terms of x, M , g~s 

• - F « {1 -^ ml = 2Mf-^±. (193) 



8B A/ (3-2i)' ' " " g s x 
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Evidently there are two combinations of eqs. (|192|) , ( |193| ) which do not depend on M 

^^^Mihg), (194) 

p = g s x(3 - 2x) = 5 u [ I . (195) 

The product gsx is positive for solutions delivering a true minimum of the effective 
potential both in the subcritical and overcritical regimes (see the discussion in section 



4). Hence x < 3/2. For the central value of L 5 = (2.2 ± 0.2) x 10~ 3 |28j one gets 

5 ~ 0.56, p~0.12. (196) 

We now eliminate gs from the mass-gap eq. ( p. 10| ) in favour of p 

(2 - p)x 2 - 5x + 3 = 0, (197) 

finding 

Clearly if p < 2 then both solutions are positive but only the second one, x_ < 1, 
delivers a true minimum. It corresponds to the subcritical regime with the positive 
gs- The first solution x + > 3/2 does not correspond to an absolute minimum. In 
the QCD case we have estimated that p ~ 0.12 and therefore the ECQM effective 
action motivated by QCD is realized in the subcritical domain. This is one of our 
main conclusions: 

x~0.92, # 5 ~0.11>0, (3>0. (199) 

Of course the fact that gs <C 1 does not by itself imply that the coupling constant 
associated to the NJL-type operators is therefore relatively weak as gs is subject to 
an arbitrary normalization. However, we have already argued when discussing the 
gap equation that the normalization used here is the natural one. Indeed that the 
theory is in a weak coupling regime can be seen in that So is close to Mq. 

Using eq. ( |192| ) one finds the central values for the bare constituent mass, M ~ 
190 MeV (i.e. S ~ 180 MeV) and consequently J ~ 0.1; that is A ~ 1.4 GeV, 
in very good agreement with the naive estimate from ChPT, A ~ 1.2 GeV. For the 
above values of B , x, M , Io one finds from (|152[) that c\ ~ — (185MeV) 3 . 

Using the relation for the heavy pion mass ( |123| ) we can estimate the asymmetry 
between scalar and pseudoscalar four-fermion interaction 

Thereby the pseudoscalar four-quark coupling, gp ~ 0.16 has a tendency to be 
stronger than the scalar one. However the pretty large errors in the determination 
of the phenomeno logical coupling constants L^,L§ makes this conclusion somewhat 
qualitative. 

Finally we apply eq.( |174| ) to fix the parameter e. According to ( |1 74|) it may have 
two values: e ~ —0.14 or e ~ —0.7. The former one is close to zero, compatible with 
a hidden left-right symmetry of the ECQM action fj\ whereas the latter one is quite 

9 The exact symmetry is realized for a scalar meson with a mass ~ 650 MeV 
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compatible (within the error bar) with the possibility to have only one type of the 
mass term, e = ±0.5 . 

After substitution into the definition ( |157| ) of the weak decay constant of the heavy 
pion Fji one arrives at the prediction 



^ _ F n ml = VT=¥ ( 2e 7 M | ^ 
' F ml 5 \g s xg A Bo / 



v/I ^ ~ ±1.0 (201) 



The constant i*n is not yet measured experimentally but has been estimated recently 
using modified sum rules PS| to correspond to \di\ = 2.2 ± 0.8. Finally, we estimate 



the scalar meson decay constant F a ~ lAF n , which quantifies the deviations from 
the predictions of the linear sigma model. 

By now we have derived the following three predictions: m a (for a given value 
of Lg, or the other way round), Fu and F a . We have so far considered the central 
values. Let us now examine the error bars. A conservative estimate 400MeV < m a < 
1200MeV is given in [11]. From fll79p it leads to 0.7 x 10" 3 < L s < 3.5 x 10" 3 , to be 
confronted with the chiral dynamics estimate L 8 = (0.9 ± 0.4) x 10~ 3 . We see that 
too light values of m a are ruled out from pion and kaon physics since the allowed 
range for L 8 leads to 700MeV < m a < 1200MeV. Conversely, one sees from ( |179| ) 



that L 8 > 0.6 x 10 3 Then, assuming compatibility of these phenomenological data, 
we derive 

L 8 = (1.0 ± 0.3) x 10" 3 , 700MeV < m a < 1200MeV. (202) 
For these ranges we find from ( |201| ) and ( p.76| ) 



\d x \ =0.6-2.5, F n = (0.4- 1.8) x lO^i^, F a = (1.2 - 2.6)F ff . (203) 

These quantities have not yet been measured. It is worth pointing out that they 
are sensitive to the sigma meson mass and therefore constitute a good laboratory for 
improving our knowledge of the scalar channel. 

We have seen that all parameters of the ECQM can be derived from meson phe- 
nomenology at low and intermediate energies and that this model is predictive. It is 
quite remarkable that the values we obtain for the different coefficients and constants 
of the effective action are in excellent agreement with our theoretical expectations. 
For instance, fo is smaller but of the same order as F ; the coupling constant gp\ is 
slightly bigger, but again of the same order as g$o, as some theoretical considerations 
suggest (see appendix B). The chiral symmetry breaking scale A comes out naturally 
in the expected range, while one of the allowed values of the parameter e is compatible 
with the left-right symmetry discussed in the text. 

9. Conclusions 

Let us summarize the key points of our approach, outline the main results and 
explore possible continuations of our work 
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1) We have seen that the original chiral quark model is not sufficient to describe 
meson physics in the GeV energy region. The structural constants L 5 , L s of the effec- 
tive chiral lagrangian also require four-quark interaction operators to be adequately 
accounted for. Moreover it can be shown (see appendix B) that a consistent QCD 
bosonization, with the help of bosonic chiral fields incorporated into the constituent 
mass operator, generically induce four-quark interaction of the NJL type. Thus the 
most notable features of the ECQM derive quite naturally from QCD. 

2) We have proposed a very general approach based on the wilsonian ideas of 
retaining only the light degrees of freedom and using symmetry and dimensional 
considerations. Thus, the ECQM should not be regarded as 'just another model' for 
low energy QCD, but rather as a genuine and correct parametrization of QCD below 
the CSB scale A in the scalar and pseudoscalar channels, much in the same way as 
the chiral lagrangian is not a model of QCD below the meson threshold, but QCD 
itself. 

3) Two novel ingredients, namely the matching to QCD, via the chiral symmetry 
restoration at high energies, and the renormalization group analysis play a crucial 
role in eliminating some of the unknown parameters and setting the accuracy of our 
approach. 

4) We have made some approximations such as neglecting quark operators of 
canonical dimension d > 6 and purely chiral operators with four derivatives of chiral 
dimension d > 4, the large- N c approach and the large-log approximation. These sys- 
tematic approximations have obviously some bearing on the accuracy of the present 
approach to the ECQM, but it is hard to go beyond them without more precise knowl- 
edge of the value of the bare couplings at the scale A, something which is at present 
impossible without solving QCD (perhaps by performing numerical simulations in the 
ECQM and matching them with QCD[[3(J). At any rate we do not really expect the 
uncertainty of the present approach to be worse, at comparable energies, than the 
one characteristic of ChPT. 

5) The approximations we have made in our action are theoretically self-consistent 
to ensure renormalizability of the low-energy action. It is remarkable that the number 
of parameters which have been retained happens to be sufficient to fit well the phe- 
nomenology of it, II and a mesons and this with quite natural values as dictated by 
dimensional and symmetry considerations. The model allows to reliably estimate the 
scalar meson masses and decay constants from other data. They come out naturally 
in a range which is not in obvious conflict with low energy phenomenology, something 
that does not happen in the NJL model. 

6) Our analysis shows that the incorporation of the chiral constituent mass Mo, as 
in the conventional CQM, leads to a relatively small value for the coupling constant 
of the effective NJL-like four fermion interactions, which come out to be relatively 
weak. 

7) It is remarkable that the scalar and pseudoscalar fermion self-interaction do not 
appear to be equal. This is because chiral symmetry is realized in low energy QCD in 
a rather subtle and not totally manifest way at the level of the effective action. It is 
worth pointing out also that the structure of soft breaking mass terms is compatible 
with the one corresponding to the L <-> R symmetry in eq. (136). 
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8) The generalization to three light quarks is in progress. One adds then two 
more 4-fermion operators to the list ( |65|) - (|68|) , but gains the physics of the scalar 
and pseudoscalar SU (3)-octet mesons. 

9) In this spirit we propose the pseudoscalar singlet state to be the natural heavy 
partner of the a meson. Thereby the resolution of the {7a(1) problem and the rj' 
meson mass will be implemented in the ECQM by means of 4-quark flavour-singlet 
pseudoscalar forces. Obviously one should go away from the large N c limit to study 
this extension. 

10) The generalization to include scalar isotriplet or spin > 1 quark currents is 
quite straightforward. The effects of such resonances can be partially absorbed in the 
structural constants of the ECQM, for instance, in 8g^. Hence their incorporation 
will modify also the latter ones (a related discussion in the case of extended NJL 
model can be found in || j26|). 
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Appendix A 

Let us justify here the use of the EOM ( |TT|) and (|23|) to eliminate certain higher- 
dimensional vertices in section 3. For simplicity we do not introduce external sources. 
The simpler case concerns the EOM of quark fields (|23|) as they are linear. The CQM 
leads to the generalized Dirac operator 

Wcqm = i$- M (U + P R + UP L ), (Al) 

acting on quark fields q(x). Since we eventually integrate out the quark fields, we 
need to use the field equation Dcqm<i{x) = inside the fermion loop to eliminate 
the unwanted operators. We shall treat all the operators with dimension > 4 as 
perturbations and thus consider the v.e.v. of products of them, performing the ap- 
propriate Wick contractions. The CQM Dirac derivative will thus act on the fermion 
propagator 

q(x)q(y) = (T(q(x)q(y))) = <ar|^ — |y>. (A2) 

P>CQM 

giving 

Dcqm q(x)q(y) = 5 (4) (x - y)I. (A3) 



The Dirac derivative ft appears in operators of d > 5 like (|T5|) , (|26|) , fl2~T|) , ([43] 
etc. The general structure is 



d xq(x)O n ftq(x) or J d xq(x) ftO n q(x) , (A4) 

where O n is a differential operator of finite order and dimension n > 1 which contains 
also the relevant matrix structure, in particular, left or right projectors if necessary. 
Thus one gets, for instance, 

J d 4 x J d*y q(x)O n ft q(x)q(y) O m q(y) = (A5) 

-Jd 4 xJ d%q(x)O n M (U + P R + UP L ) q(x)g(y) O m q(y) + J d A x q(x)O n O m q(x) 

The last operator in (A5) has dimensionality d = n + m + 3 > 5. Thus the application 
of the EOM generates operators of at least the same dimension without the Dirac 
derivative (or, in general, with one Dirac derivative less). The result of the application 
of the EOM is the removal of all Dirac derivatives acting on quark fields at the expense 
of modifying the coupling constants in operators of d > 5 without derivatives. In 
particular, the d = 5 operators are mapped into a set of d > 5 operators, the d = 6 
vertices are mapped into a set of d > 6 operators, etc. 

One can improve this reduction when dealing with the ECQM effective action 
(ff4l). Indeed, the latter one includes all vertices of d < 4 and therefore the higher 
dimensional operators start now from d = 5. It corresponds to m > 2 in (A5). Thus 
in the ECQM, with the help of the EOM, the d = 5 vertices are mapped onto a set 
of d > 6 vertices etc. 

In an analogous way one can show that the use of the EOM for chiral fields ( |TTD 
leads also to the modification of coupling constants in higher dimensional vertices 
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without an insertion of the operator d^. The proof is now complicated by the 
nonlinear structure of the EOM. Nevertheless one can expand the chiral variables 
U in powers of 7r and analyze the pertinent EOM reduction in ChPT. The ChPT 
counting rules assign dimension zero to the pion field, and ~ d 2 n thus carries 
chiral dimension d = 2. The equations of motion can clearly be applied to external 
pion fields without further ado. As for those operators contributing to internal pion 
loops they will give rise, after the appropriate contractions, to the pion propagator 

(T(n(x)rr(y))) = (x\±\y). (A6) 

The application of the Klein-Gordon operator d 2 arising from the EOM to the above 
Green function gives 

d 2 (T(n(x)n(y))) = 5^(x-y)l. (A7) 

Since all operators of d < 3 are included in the effective lagrangian of the CQM the 
set of higher dimensional operators treated as perturbations consists of those ones 
with d > 4. Hence after the EOM reduction in two vertices connected by one pion 
propagator one expects that additional local vertices of g? > 4 + 4 — 2 = 6 are induced 
similarly to (A3). Evidently these d = 6 operators have also chiral invariant structures 
and therefore represent particular combinations of the complete set of d = 6 operators 
classified in section 3. If some of them still contain insertions of the operators 
and/or ft acting directly on quark fields, repeated use of the EOM produces additional 
operators of d > 6 without any such insertions. 

Appendix B 

Below the chiral symmetry breaking scale A, the (in the chiral limit) massless 
Goldstone bosons, associated to the local chiral phase of the quarks, appear. To 
avoid overcounting we must extract the chiral phase from the QCD partition function. 
This naturally modifies the QCD action. Let us see how and that, in particular, the 
separation of chiral fields in QCD is in general accompanied with the pertinent four- 
quark interaction. We separate the local chiral phase by exploiting the following 
functional identity (for a related approach see [31]) 



1 = A (VU l[S(q R Uq L (x) + q L U^q R {x)) , (Bl) 

where U are SU(2) fields and A is the required Faddeev-Popov determinant. 

Now we proceed to calculate the determinant and use the conventional mapping 
of chiral SU{2) fields onto S 3 

U = n°I + ir a n a , a = 1,2,3 {n k f = (n ) 2 + n a n a = 1, (B2) 

with the notation k = 0, 1, 2, 3. Then the properly normalized integral over chiral 
fields is evaluated as follows 

/ VU\{5{q R Uq L {x) + q L U^q R {x)) = J H^ls((n k (x)) 2 - 1) 5(J k (x)n k (x)) 
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= I[(Jk(x)J k (x)r 1/2 , (53) 

X 

where C is a normalization constant and the quark currents are 

J (a:) = qq(x); 3 a (x) = iq^r a q(x). (BA) 
Therefore the determinant supporting the identity (Bl) reads 

A = J] ((qq(x)f - {qi,r a q{x)f) 1 ' 2 . (2?5) 

X 

This four-quark operator can be exponentiated with the help of auxiliary fields. The 
S function in (Bl) can be exponentiated as well by using the functional Fourier rep- 
resentation. All together allows us to write 

1 = J VUV[iV(f)VcVc exp J d^x (ifx(x)(q R Uq L (x) +q L U ] q R (x)) 

-(0 2 (x) + cc(x))((qq(x)) 2 - (qi 5 T a q(x)) 2 )), (BQ) 

where U(x), /jl(x), <f>(x) are bosonic fields and c(x), c(x) are scalar fermionic ones (ghost 
fields). 

After insertion of this identity in the QCD generating functional one can argue 
that the quantum average over the gluon vacuum should generically induce (since they 
are not protected by any symmetry) non-zero v.e.v. for the auxiliary fields (/i) = Mq 
and (4> 2 (x) + cc(x)) = g/4iV c A 2 . Should all the above v.e.v. be equal to zero we 
would conclude that the chiral symmetry is not broken and the U field would entirely 
decouple. Although we cannot at present demonstrate the actual mechanism, we 
know that this is not the case in low energy QCD and from that point on the analysis 
presented in this paper takes over. 

Appendix C 

Here we shall compare our results with those obtained in the pure CQM after 
inclusion of gluonic corrections and see how they fit together. In the values for Fq 
and the low energy chiral constants Lj were derived considering both the contribu- 
tion which is radiatively induced by the fermion loop with constituent mass M and 
additional gluon contributions expressed in terms of v.e.v. of gauge invariant gluon 
condensates. The leading gauge invariant operator is naturally (a s GG) and this was 
the only one explicitly considered in . 

As discussed in the text, the chiral constants L 1; L 2 , L 3 , L 9 and L w , which in the 
CQM get a contribution from the constituent quark loop 

N 

8Lt = 4L 2 = -2L 3 = L 9 = -2L W = — (CI) 

48tt z 

are exactly reproduced by the ECQM. The gluonic corrections obtained in would 
correspond to what in this paper we call bare contributions (and which we have 
explicitly considered here only for F ). For the present discussion we take g& = 1, 
just as the authors of Q did. We have obtained the following expression for F 

F 2 = N c M 2 h + fl (C2) 
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With the same notations the following result was obtained in 

F 2 = N c Mlh + ^mI^sGG), (C3) 

and the identification of /o is obvious. Notice that the scale \x which was undetermined 
in [[J is to be identified with the CSB one A. We know this because we have good 
control on the renormalization group behaviour in the low energy effective theory. 
We see clearly that the bare coefficients contain all the gluonic contributions which 
accumulate all gluonic modes both above and below A which cannot be expressed in 
terms of the constituent mass Mq. 

As to the chiral constants L 5 and L 8 , using the large log approximation, we read 



from 



and 



J , (C5) 



N c Mo 



16B 

both of which can be made to agree with our results in the gso — gpi — limit for 
5fo = 0. This limit is quite straightforward to take in L 5 (see (|189|) ), but it is worth 
recalling that in the present approach L 5 is generated through a exchange and then 
the corresponding amplitude has the correct high energy limit. The situation for L 8 
is a bit more subtle. The limit in ( |179| ) cannot be taken naively in the final expression 



for Lg which vanishes, but rather we must replace Z a and Zji by their expressions in 
the first term in ( |179| ) and then expand in powers of gso = 9pi- We then recover (C5). 
The reason for this subtlety must be eventually traced back to the CSR constraints. 
On the other hand, in [TJJ, and also using the leading log approximation, one gets 

C q = N c MqI (C6) 

This shows that within this approximation M , and consequently L 5 , must be neg- 
ative, in contradiction with experiment. In this difficulty is circumvented due to 
the presence of (unfortunately non-universal) finite parts which contribute to L 5 with 
the opposite sign and eventually allow getting around this sign problem, although 
the final fit for L5, L§ is not brilliant. We have seen in the main text how in the 
ECQM the renormalization group implies, on the one hand, the presence of some 
logarithmically enhanced contact terms which make possible for C q and L5 to have 
opposite signs without having to appeal to finite non-universal constants. 
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